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Abstract 

We first show that the tihngs of a do main D form a lattice (using 
the same kind of arguments as in |Rem99|) which we then undertake to 
decompose and generate without any redundance. To this end, we study 
extensively the relatively simple case of hexagons and their deformations. 
We show that general domains can be broken up into hexagon-like parts. 
Finally we give an algorithm to generate exactly once every element in 
the lattice of the tilings of a general domain "D. 
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1 Introduction 

Tilings is an age-old topic for specialists and amateurs alike. Lozenge tilings in 
particular have intrigued generations of curious people because they can easily 
be seen as piles of cubes. 

In the past few years, mathematics, theoretical physics and computer science 
have started shedding a new light. We will build on these results to provide an 
answer to a very natural question: given a tileable domain and sufficiently 
many lozenges tiles, how can one generate all the tilings of the domain without 
repeating twice the same tiling? 



Conway and Lagarias introduced in |CL90| a new, powerful tool to study 
tilings-related topics: tiling groups, which give a necessary condition for a do- 
main to be tileable, and provide and important bijection in the case of lozenges. 



Thurston went a step further in |Thu90| and showed by a constructive algorithm 
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(a) The two admissible 
regular triangles 



(b) There are only three ad- 
missible lozenges 



Figm'e 1: Triangles and lozenges 



that one can quickly decide whether a domain is tileable; to this end he uses 
height functions (see Section |[). 

Thurston also hinted that the set of the tilings of a domain, partially ordered 
with the height functions, should have the structure of a lattice. This was proved 
in |Rem99 by Remila in 1999 in the case of dominoes. We adapt his proof to 
the case of lozenges in Section |[ 

We then proceed to the really new material: after a rather extensive study 
of the case of hexagonal and hexagonal-like domains (see Section |^), we use 
a geometrical point of view (justified by a bijection between Conway and La- 
garias' lozenge group and Z^) and the identification of meaningful hexagonal- like 
sub-domains (see Section ^ to exhibit in the general case a maximal chain of 
intervals in the lattice of the tilings (see Section |^). This chain is a natural 
extension to Thurston's minimal and maximal tilings. We finally introduce new 
minimal tilings which allow us to achieve our goal: (uniquely) generating all the 
elements of the lattice formed by the tilings (see Section ^. 

We believe that the tools introduced in this paper, notably the seeds and 
their ranges, should prove fruitful in tackling related problems (see Section p|). 



2 Basic tools and definitions 

In this section, we present the definitions of classical objects which will be used 
in this paper, in an attempt to make it reasonably self-contained. All the new 
objects are defined later in the paper, at the moment when they are needed. 



2.1 Tiling with lozenges 

First of all, let us define what we mean by a tiling. One needs two regular 
triangles (see Figure (a)). The whole plane can be covered with a repetition 
of these figures (see Figure ^ (a)), which gives rise to the triangular grid. A 
domain is a finite union of triangles in the grid. It is simply connected if it is 
connected and its complement in the plane is connected. A polygon is a simply 
connected domain. 

We now define three tiles by gluing together two regular triangles; let us call 
them lozenges (see Figure |] (b)). A tiling of I? is a set of lozenges included in 
V with pairwise disjoint interiors such that the union of the lozenges is V itself. 
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(a) The plane can be covered (b) A general domain T) lim- 

with regular triangles ited by a closed path V 



Figure 2: A closed path in the triangular grid 




A domain is tileable if it admits at least one tiling. Its boundary or contour is 
the set of its edges that belong to exactly one of its triangles. 

Any finite-length closed path V whose edges belong to triangles in the trian- 
gular grid can be viewed as the boundary path of a domain V (see Figure || (b)). 

2.2 Tiling groups 

Figure H (b) suggests a connection between tilings by lozenges and certain piles of 
cubes, which we now attempt to clarify. Our tool here is Conway and Lagarias' 
tiling groups, a formal description of which can be found in [CL90, rhu90 , so 
we will restrict ourselves to an intuitive (but rigourous) approach. The idea is 
to simplify contour words so that those of tileable domains are equivalent to the 
empty word. 

First, label the sides of triangles with letters; in our case, the set {a, b, c} is 
sufficient (see Figure ^ (a)). Our purpose is to use words to partially describe 
tilings. Let T be a tiling of a polygon and let P be a directed path (in the 
triangular grid) whose edges belong to lozenges in T; that is, P never cuts 
the interior of a lozenge. Such a path will be called T-valid or valid for T. 
Given a starting point, P is completely encoded by a word x on the alphabet 
{a, b, c, a~^, b~^, c-^}. 

We want P to be rather straightforward, so that following an edge and then 
following it in the opposite direction should not change x. We thus impose 
aa~^ = bb^^ = cc^^ = e (the empty word) and make all possible simplifications 
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(a) 



(b) 

Figure 4: Labelling the sides 



in X, so that P is now encoded by a word in the free group F{{a, b, c}). Although 
this is not a pre-requisite for the rest of this paper, the reader unfamiliar with 



free groups can look up page 257 in ||Fra94| for instance 



If P' is another T-valid path, encoded by a word y in F{{a, b, c}), with the 
same starting and ending points as P, then P and P' define a polygon V which 
is tiled by lozenges, and therefore tileable. To reflect this fact, we set the contour 
word xy~^ of 2? to be the empty word, e. In particular, the contour words of each 
of our three lozenges should be set to e, that is (see Figure ^ (b)) bcb~^c^^ — e, 
aba~^b~^ = e and cac~^a~^ = e by reading counterclockwise (starting from a 
different vertex only changes the word by a circular permutation). Note that 
this process is equivalent to removing loops in P, or stating that two T-valid 
paths having the same starting and ending points should be described by the 



same word, as we will see in Corollary 2.4 



The words we consider now belong to the group 

L ^ (a, 5, c I bcb^^c^^ = aba^^b^^ = cac^^a^^ = e) 
which is Conway and Lagarias' lozenge group. 



Definition 2.1 Let T be a tiling of a polygon and let P = {vq, vi,. . . , Vp) be 
a valid directed path in T. Each edge (vi,ViJ^i) can be labelled by an element 
of the alphabet {a, 5, c, a^^, 6^^, c^^}; the label of P is the word w obtained by 
the concatenation of the labels of its edges in the order (wq, vi), . . . ,(wp_i, Vp). 
The free label a{'w) of P is the representative element of w in the free group 
-F({a,5, c}). The L-label £{w) of P is the representative element of w in the 
lozenge group L. 

For the sake of simplicity, the free label and the L-label of P will also be 
written a{P) and £{P)- 

The following lemma gives a necessary condition for a domain to be tileable: 

Proposition 2.2 ( [|CL90| ) If a polygon is tileable, then the L-label of its 
boundary path is trivial. 

Proof We make the proof by induction on the surface. (The reader familiar 
with De Bruijn's worms will readily find a direct proof.) It is enough to prove 
the result for elementary cycles. Besides, the result holds for single lozenges. 
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Figure 5: Projection of a cube onto the plane x + y + z = 



Let P be a valid path in a tiling T of the polygon V, such that P is distinct 
from the boundary path B oiV . P cuts V into two tileable polygons Vi and 1^2- 
There exist w and w' such that the boundary paths of Vi and 7^2 ai'e 'wa{P) 
and q;(P)~^z«', both of which are trivial by the induction hypothesis. Therefore 
i{B) = £{w ■ w') = £{wa{P) ■ a{P)-^w') = £{wa{P)) ■ £{a{P)-^w') =e. □ 

Corollary 2.3 Let T be a tiling of a polygon. The L-label of any closed valid 
path in T is trivial. 

Proof Indeed, a valid path delimits a tileable polygon of which it is the bound- 
ary path. □ 

Corollary 2.4 Let T be a tiling of a polygon. Two valid paths in T having the 
same starting and ending points have the same L-label. 

Proof Indeed, if w and w' are the labels of these paths, w^^ ■ w' is the label 
of a valid closed path in P. □ 

The relations that appear in the definition of L can be rewritten as ab ~ ha, 
ac — ca and be — cb, so that L has three generators that commute with each 
other and is therefore isomorphic to Z'^. How can we interpret this nice result? 

Let P be a tiling of a polygon V and let w be a vertex of P. A valid path 
in P can be associated with a unique word in L, which in turn corresponds to 
a unique path in the 1-skeleton of a cubical tesselation P of space. There is 
therefore a one-to-one correspondence between vertices of V (resp. segments 
in P) and vertices of Z"^ (resp. segments in P). Following an edge in P is 
analogous to following one of Z'^'s generators. This bijection allows us to lift 
each edge of P to an edge in P, each lozenge to a square in Z^^, so that P is 
equivalent to a collection of squares in P. This squares may or may not define 
the visible parts of cubes, depending on P. The squares can be projected along 
the direction (1, 1, 1) in to give back P (see Figures ^ and ^ (b)). 



5 




(a) Lemma 2.2 does not give 
a sufficient condition 



(b) Height functions take 
integer values 



Figure 6: The need for height functions 



2.3 Height functions 

The contour word of a polygon may well have a trivial image in the lozenge 
group even though the polygon is not tileable. In other words, the condition 
stated in Proposition ^.2| is not sufficient. Consider for instance Figure (a): 
the polygon is clearly non-tileable, but its contour word hccab~^c~^c~^a has 
a trivial image in L because letters commute. More generally, any closed path 
in 1? corresponds to a word in which the number of a's (resp. 6, c) is equal 
to the number of a~i's (resp. c~i), so that the image of this word in L 

is trivial; it would be quite surprising if the projection onto the plane of any 
closed path in 1? was tileable. More complicated examples can be exhibited, 



for instance using Fournier's obstructions (see Fou96 1 ) . Instead of examining 
each case, let us develop from this simple example a general idea. 

Let a, b and c correspond to (1,0,0), (0, 1,0) and (0,0, 1): we can now see 
Figure || (a) as a closed path in 1? (it looks like a cyclohexane molecule in chair 
conformation). If the domain was tileable, one could view a tiling as a collection 
of squares and cubes in I?. Since the distance between A and B is only one 
edge in the triangular grid, these points would be linked by either the side or 
the diagonal of a square in Z^, so that their distance in would be 1 or ^/2; 
but if A has coordinates (0,0,0) then B has coordinates (1, 1,2) since starting 
from A we follow b once, c twice and a once. 

The distance (in 1?) between A and B can easily be made greater: remark 
that the figure looks like a butterfly, and make the wings bigger without moving 
A or B. This adds 3 to the distance between A and B at each step. 

We see that the distance between two points seems important. Two points 
at distance 1 in the triangular grid should not be distant ones in 1? if the 
domain is tileable. What really matters in our butterfly example though is 
the distance between A and B along the (1, 1, 1) axis: following an edge in a 
tiling is like following one of Z'^'s generators, each of which yields the same 
height increase, say, i, along the (1,1,1) axis so that every point in Z^, once 
orthogonally projected onto the (1, 1, 1) axis, is at a distance to the origin that 
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is a multiple of z. To each point p of one can therefore give a value (multiple 
of i) which corresponds to the distance between a fixed origin point and the 
orthogonal projection of p onto the (1, 1, 1) axis. For convenience, we will forget 
the geometrical interpretation, place the origin at an arbitrary vertex and set i = 
1 (see Figure |^ (b)). 

We now proceed to properly define height functions, mainly following Thurs- 



ton (see |Thu9C|) while preserving the algebraic point of view. 



Definition 2.5 We call evaluation function the morphism (p from the lozenge 
group (L, •) with generators a, b, c, to the group (Z, +) of integers such that 

ip{a) = (p{b) = ip{c) = 1. 

Note that this implies (p{e) — and tp{w ■ w') = ffiiw) + f{w') for any w, w' 
in L, whence Lp{a~^) = ip{b^^) = ^p{c^^) = —1. 

Definition 2.6 Let T be a tiling of a polygon V and let w be a vertex on the 
boundary path of V . The height function induced by T and v is the function 
that maps each vertex a; of 7^ to the image by the evaluation function of any 
valid path from v to x. 



The correctness of this definition stems from Corollary 2.4 



Lemma 2.7 On the boundary path of V, the heights of the vertices do not 
depend on T. 

Proof Indeed, if x is such a vertex, there exists a path from u to x that lies 
on the boundary path oiV , so the result follow by induction. □ 

Note that only the vertices on the boundary path have fixed heights; those 
of inner vertices do depend on T (see for instance Figure |^ (a)). 

Since changing the reference vertex v only changes the height function by a 
constant, we will often refer to a height function without mentionning v. Thus 
we will write "the height function associated with the tiling" . 

2.4 Thurston's algorithm 



We have seen (see Section 2.3) that Proposition p.2| does not give a sufficient 



condition for the tileability of a polygon. Thurston's height functions provides 



a constructive algorithm, outlined in |Thu9C|, to determine whether a polygon 



can be tiled, and exhibit a tiling if it can be done. 

We will build the minimal tiling of the polygon, in a sense that will be made 
clear in Section ^. There is a natural (partial) order on the height functions as 
the tiling changes and this is interpreted as an order on the tilings. 



The following algorithm, defined by Thursto n in | Thu9C ] , is based on this 
simple result (see Proposition 3/7 and Definition 3^ in Section || to know more 
about flips): 
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Lemma 2.8 Let v be a vertex of a polygon V such that the height function 
associated with the minimal tiling of V is maximal on v. This vertex cannot 
belong to the interior of V, otherwise it could be flipped down: therefore v lies 
on the boundary path of V . 

Algorithm 2.9 

• Input: A polygon V . 

• Output: The minimal tiling of V if the polygon is tileable, untileability 
otherwise. 

• Initialization: Initialize the list L to 0. 

• Step 1: If 7-" is a single point, return L. 

• Step 2: Compute the height function on the boundary path of V. If a 
vertex is given two different heights, return untileability. 

• Step 3: Let w be a vertex of V of maximal height. There exist u and 
w on the boundary path B of V such that {u,v) and (v,w) are segments 
of B. Place a lozenge i on V so that u, v and w belong to it. Add £ and 
its position to L. Update V to V \ £. Go back to step 1. 



For the proof of the algorithm, the reader can refer to | Thu90| . Its complexity 



is linear in the number of triangles in V. From this algorithm one can easily 
deduce an algorithm to build the maximal tiling. 

2.5 Lattices 



We will see in Section B.2 that the height functions associated with the tilings 
of a polygon define a partial order on this tilings; this order has the interesting 
property of being a lattice, which we now define. See Figure ^ for an example 
of a graphical representation of a lattice. 

Definition 2.10 (Lattice) A set S partially ordered by a relation is a lattice 
if for any a and b in S there exist i and s in S* such that 

• i ^ a ^ s and i ^ b ^ s; 

• if X ^ a and x ^ b then x ^ i for any x € S] 

• if a ^ y and b ^ y then s ^ y for any y G S*; 

Such elements i and s are called the infimum and the supremum of a and fe, 
noted inf(a, b) and sup(a, b) (or a A b and a V & for short). If a and b happen 
to be comparable (e.g. a ^ b) then i and s are called the minimum and the 
maximum of a and b, noted min(a, b) and max(a, b) (equal respectively to a and 
6 if a b). 
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(a) The Ferrers diagram of the (b) A plane partition 

partition (4,4,2,1) of 27. 

Figure 7: Ferrers diagrams and plane partitions 



If 5* is a finite lattice (such as the tiling lattices considered in this paper), 
then it admits one maximal and one minimal element. 

A lattice {S, ^) is distributive if for any a, b and c in 5 one has 



a A (6 V c) = (a A 6) V (a A c) 
aV (6Ac) = (a V&) A (a Vc). 



Definition 2.11 (Interval) An interval [a;b] (with a 6) in a lattice (C,^) 
is the set of all a; G £ such that a ^ x ^ b. 

Note that a and b need be comparable elements. 

Proposition 2.12 Let (£, =4) be a lattice and T an interval in £. The order =^ 
defines on X a structure of lattice. 

2.6 Partitions, Ferrers diagrams and plane partitions 



These definitions will be needed in Section 4.2. See also Figure 



Definition 2.13 A partition of an integer n is a non-increasing list of positive 
integers with sum equal to n. 



Definition 2.14 A Ferrers diagram is a geometrical representation of a par- 
tition: given a partition (Ai,... ,Ap), the corresponding Ferrers diagram is a 
collection of p consecutive left-justified rows respectively containing Ai, ... , Ap 
consecutive squares. 



Definition 2.15 A plane partition of an integer n is a filling of a Ferrers di- 
agram with integers (of sum n), such that the values along the rows and the 
columns is non-increasing. 
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(a) A hexagon can be 
tiled in two ways 



(b) Height functions for the 
three fozenges 



Figure 8: Our building blocks 



3 Flips and lattices 

In this section, we first show that there is a one-to-one correspondence between 
tilings and height functions (up to constant) for a given polygon V , which allows 
us to show that the set of the tilings of V under the partial order defined by 
the height functions has the structure of a lattice. We then consider a local 
operation (called a flip) on a tiling and prove that it has a straightforward 
translation in terms of the underlying height function. 

A basic remark will be the leading idea for most of this section: a hexagon 
can be tiled in exactly two ways with lozenges (see Figure || (a)). Switching 
between these two tilings is called flipping. This operation allows us to deduce 
new tilings from a known one. The really challenging idea is that this allows us 
in fact to navigate between all the tilings of T). All the results and proofs in 



this section are adapted to the case of lozenges from |Rem9£], which deals with 
dominoes. 

3.1 Prom height functions to tihngs 

We investigate in this section the link between height functions and tilings. An 
edge will be said positively directed if it is labelled by a, b or c. Recall that an 
edge in the triangular grid is valid for a tiling T if and only if it does not cut 
the interior of a lozenge in T. 



We have seen (see Definition 2.6) how a height function h can be defined, 
starting from a particular tiling T of a polygon V. Along any valid positively 
directed edge the height changes by 1; along any invalid positively directed edge, 
the height changes by —2 (see Figures ^ (b) and || (b)). Therefore h is merely 
an encoding of T , and the latter can be reconstructed from the former. 

Lemma 3.1 There is a one-to-one correspondence between the tilings of a poly- 
gon V and the associated height functions (up to constant). In other words, if 
T and T are two tilings of V and if hxiv) = hr'iv) for any t; e P then T = T' . 

Proof Assume V and a height function Ht given, the latter corresponding 
to an unknown tiling T oi V. To rebuild T from hx, it suffices to draw all 
the triangles inside V and to erase those edges whose endpoints have a height 
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difference of 2 in absolute value. These edges lie on the interior of lozenges and 
completely characterize them so we are done. □ 



Proposition 3.2 Let T and T' be two tilings of a polygon V. For each vertex 
V oiV , hxiv) — hT'{v) is a multiple of 3. 

Proof We make the proof by induction on the vertices. The proposition is 
true on the bound ary path of V since all height functions take the same values 



there (see Lemma 2/7). Assume that the proposition holds for v and let v' be 
any of its neighbours in the triangular grid. Let T be either T or T'. If {v,v') 
is a valid positively directed edge in T then hrW) — h-r{v) — 1; otherwise, 
hrW) — hr{v) = —2 and thus h-r{v') — hT{v) takes the same value modulo 3 
whether {v,v') is valid or not. Consequently, if hxiv) — hT'{v) is a multiple of 
3, then so is hriv') — hT'{v'). □ 

An intuitive way to consider this proposition is to think in terms of cubes: 
adding 3 to the value of a height function on a vertex v is equivalent to adding 
exactly one cube (see Figure ^ (a)). 



3.2 The lattice structure 

In this section, we prove that the set of the tilings of a polygon can be endowed 
with a lattice structure. To this end, we define a partial order between tilings 
by using their height functions: 

Definition 3.3 Let hr and hx' be two height functions associated with the 
tilings T and T' of a polygon V such that hr and hx' take the same values 
on the boundary path of V. hx is less than hx' (and we note hx ^ hx') if 
hx{v) ^ hx'{v) for any vertex v dV. 

The order between height functions allow us to endow the set of the tilings 
of V with a natural order: T =<; T' if and only if hx ^ hx' ■ 



Proposition 3.4 Let T and T' be two tilings of polygon P and let hx and 
hx' be their height functions. The functions /imin = niin (/i^, hx') and /imax = 
max {hx, hx') are themselves height functions. 

Proof We prove the result only for /imin since the other case is analogous. We 
make the proof by induction. Since hx and hx' take the same values on the 



boundary path B of V (see Lemma 2.7), /imin(w) = hx{v) = hx'{v) for every 



V & B. Let {v,v') be any positively directed edge between vertices of V in the 
triangular grid. We claim that hyain{v') — ft.min('f) equals either 1 or —2. 

• If hx{v) ~ hx'{v) then /imin(v') — ft-min(w) cquals cither hx(y') — hx{v) or 
hx'{v') — hx'{v) and must therefore be equal to either 1 or —2. 
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A 



Figure 9: An example of lattice 



• We can now assume without loss of generality that hriv) < /it' (v), so that 
hmin{v) = hriv). hxiv') can only be hT{v)+l or hT{v)—2, and Ht' (v) is at 
least /iT(f) + 3 by Proposition |3.2| . Moreover hx'iv') can only be hT'{v) + l 
or hT'{v) — 2, and therefore hT>{v') is at least hT{v) + 1. Thus h^i^W) = 
hxiv'), from which we derive hiain{v') — hminiv) = hriv') — hriv) which 
must be either 1 or —2. 

We have shown that /imin increases by either 1 or —2 along any positively 
directed edge in the triangular grid. Consequently, the set of heights (for /imin) 
modulo 3 of the vertices of any triangle in the triangular grid must exactly be 
{0, 1, 2}. Moreover the height difference along an edge does not depend on which 
of the two neighbouring triangles was chosen. Erase all edges whose endpoints 
have a difference of heights of —2. What is left is a tiling of V with lozenges; 
that is, /imin is a height function. □ 

Corollary 3.5 The order =^ induces a structure of distributive lattice on the 
set of the tilings of V. 

Proof If Ti and T2 are two tilings of a domain T) then the height functions 
hniin{Ti,T2) and /imax(?i, T2) are clearly th e infi mum and supremum of the 
height functions hr^ and /1T2 (see Definition 2.1C| ), and since height functions 



encode tilings (see Proposition |3.2| ) we have defined the infimum and supremum 
of Ti and T2. 

To prove distributivity, it suffices to consider each vertex of V, which brings 
us back to checking the relation for integers and the usual min and max func- 
tions. □ 

Figure |^ provides an example of the graphical representation of a lattice 
of tilings. Intuitively, one obtains the infimum of two tilings Ti and T2 by 
selecting only the cubes which appear in both, and their supremum by selecting 
the cubes which appear in either one. In other words, inf(Ti, T2) is encoded by 
/imin(ri,T2) and sup(Ti,r2) is encoded by /imax(T'i, 12). 
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3.3 Flips 

In this section, we introduce an elementary operation classically called a flip; 
we prove that the set of the tilings of a polygon is connected by flips. 

Definition 3.6 Let T be a tiling of a polygon V. A local maximum (resp. 
minimum) of the height function Ht is a vertex v in the interior of V such that 
hxiv) ^ hriv') (resp. hriv) ^ hriv')) for any v' neighbour of v. 



Proposition 3.7 Let T be a tiling of a polygon P. A vertex v in the interior 
of is a local extremum of the height function Ht if and only if it is the center 
of a hexagon tiled with three lozenges. 

Proof Let v' and v" be two neighbours of v so that {v,v',v") is a triangle. 
If {v',v) and {v,v") were both valid edges in T, then h{v) would be less than 
h{v') and more than h{v") (or the converse), so it would not be an extremum. 
Moreover (w',w) and {v,v") cannot both be invalid edges, therefore exactly one 
of them is. Consequently, the hexagon around v is tiled with exactly three 
lozenges. 

The converse part of the proof is obvious (see Figure || (a)). □ 

Definition 3.8 (Flip) A flip is the operation by which one switches from one 
tiling of a hexagon to the other. An up-flip increases the height while a down-flip 
decreases it. 

We now prove that flips allow us to reach any tiling of V from any other 
tiling of P. To this end, we need one more definition: 

Definition 3.9 If T and T' are two tilings of a polygon V, then the distance 
between them is 



A(T,r') = ^|/iT(t')-/^T'(^^)l 



This function is indeed a distance since it is symmetrical, satisfies separation 
(A(r,T') = if and only if T = T' by Proposition 1^) and the triangular 
inequality (since | • | does). 

Proposition 3.10 Let T and T' be two tilings of a polygon V. If T =<; T' then 
there exists a sequence (To — T,Ti, . . . ,Tn ^ T') of tilings of T) such that Tp+i 
is deduced from Tp by a single up-flip, ^ p ^ n — 1. 

Proof Assume T ~< T', let m — min{/iT(w) \vGV and hxiv) < /iT'(f)} and 
let w be a vertex of V such that hxiv) = m. The vertex v cannot belong to the 
boundary path of V by Lemma pj| and hT'{v) must be at least hriv) + 3 by 



Proposition 3.2. Let v' be a neighbour of v such that the edge (v, v') is valid for 
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/iT- hT'{v') is at least hr'iv) — 2 ^ /it(w) + 1 > hriv). Because m is minimal, 
the edge {v,v') must be directed from v to v' . Since this is true for any edge 
anchored in v that is valid in /i-r, we conclude that u is a local minimum of Ht- 
By an up-flip on v, we send T to a tiling Ti such that T ~< Ti =4 T'. Besides, 
/iTi differs from Ht only on the vertex v, so that A(T, Ti) = 3. Using induction, 
we thus build a sequence of ever greater tilings. Since the distance between two 
tilings must be a multiple of 3 (see Proposition |3.2D , there comes a tiling T„ 
such that A(T„, T') — and T„ ^ T', and since height functions encode tilings 



(see Lemma |3j), T„ = T'. □ 



Theorem 3.11 The set of the tilings of V is connex: any tiling T oiV can be 
reached from any other tiling T' of V by using only flips. 

Proof Indeed, any tiling of V can be reached from the minimal one by using 
only up-flips (see Proposition and the minimal tiling of V can be reached 



from any other tiling by using only down-flips. □ 



4 Hexagons and pseudo-hexagons 

The ultimate goal of this paper is to provide an algorithm for uniquely generating 
all the elements of the lattice of the tilings of a polygon P. As we will see in 
Section ||, some hexagonal-like subsets of P play a important role. We thus 
start with the rather simple case of V being a hexagon or a pseudohexagon (see 
definitions below). 



4.1 Piles of cubes 

It is well known in tilings lore that the tilings of a hexagon of side n are bijectively 
related to some piles of c ube s. Using C onway and Lagarias' lozenge group, we 
have explicited in Section pTs] (following [[rhu90 |) the algebraic translation of the 



geometrical intuition. The bijection between the lozenge group and Z'^ implies 
that it is equivalent to consider lozenges in a tiling or 2-cells (squares) in the 
cubical tesselation of Z"^. Until now, we have been interested in manipulating 
tilings and watching the translation in the figure. For instance, rearranging the 
three tiles of a hexagon of side 1 is equivalent to raising the height of exactly one 
point by 3: this is an up-flip. Now we look at the bijection the other way: we 
manipulate squares in Z^ to produce new informations on tilings. For instance, 
it is easily seen from Figure ^ that while using only flips suffices to generate all 
the tilings of a polygon, it also generates multiple times the same tiling. 

Not all collections of squares in Z'^ correspond, once projected onto the 



plane, to tilings in the triangular grid: see for instance Figure 11 (a). So we 
need conditions on the squares. In fact, we will build our reasoning not on 
squares but on cubes, because the elementary operation on a tiling, the flip, 
has a natural, visual and intuitive translation in terms of cubes (see Figure 



Moreover, fracture hues (see Section 5.2) can be used to get rid of squares that 
do not correspond to cubes. 
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An up-flip is equivalent to adding a cube 

(a) Neither a tiling nor a com- (b) A pseudo-hexagon delimited by 

pact pile of cubes three Ferrers diagrams 

Figure 11: Compact piles and pseudo- hexagons 



The minimal tiling of a polygon has no local maximum in the interior of the 
polygon, otherwise this maximum could be flipped down. In other words, there 
is no cube associated to such a tiling. It merely provides squares on which one 
can put cubes. In the case of a hexagon, the association of the squares can be 
seen as base planes in (2+)^ so that any cube can be encoded by the (integer) 
coordinates of its lowest corner. 

Cubes can be piled, but if we want to preserve a tiling after projection, we 
must pile them in a way that corresponds to flips. Such piles will be called 
compact. 



Figure 10: 



Definition 4.1 (Compact pile) A pile of cubes is compact if for any cube of 
the pile at (io,jo,ko) there are cubes at («, jo^^o), {io,jiko) and {io,jo,k) with 
i, j and k ranging in [0, io], [0, jo] and [0, fco]. 

The next lemma easily stems from this definition: 

Lemma 4.2 A pile of cubes is compact if and only if any section of the pile by 
a plane orthogonal to one axis {Ox, Oy or Oz) yields a Ferrers diagram. 

It will not be enough to consider only "perfect" hexagons of size n x n x n. 
The hexagon-like sub-polygons that appear in tilings require a somewhat more 



general approach (see Section 5.5) 



Definition 4.3 (Pseudo-hexagon) A pseudo-hexagon is the domain obtained 
by starting from a compact pile of cubes and performing all the possible down- 
flips. 
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(a) Tilings are related one-to- 
one to plane partitions 




(b) A limited partition 
Figure 12: Plane partitions and limited partitions 



Note that by Lemma 4.2 a pseudo-hexagon can also be seen as a figure 
delimited by three Ferrers diagrams fitting neatly (see Figure |ll| (b)). It is also 
tileable by construction. 

Proposition 4.4 A pile of cubes in a peudo-hexagon corresponds to a tiling of 
the pseudo-hexagon if and only if it is compact. 

Proof A pile of cubes corresponds to a tiling of a domain V if and only if it 
can be generated (starting from the minimal tiling) by using only up-flips (see 
Corollary |3.5| ). An up-flip adds a cube that lies on squares that belong either to 
the minimal tiling or to an already added cube. The result follows by induction 



on the number of cubes and Lemma 4.2. □ 



4.2 An algorithm to generate the tihngs of a 
pseudo-hexagon 



Proposition 4.4 means that generating the tilings of a pseudo-hexagon is equiva- 



lent to generating compact arrangements of cubes, which is quite easier. Indeed, 



the latter are r elated one-to-one with plane partitions (see Definition |2.15| , Fig- 



ure y (a) and |Els84|) 



Definition and notation 4.5 The parts of the partition p associated with a 
Ferrers diagram F are noted p[j], and F[j] denotes a collection of p[j] linearly 
adjacent squares of 1?. A partition on F[j] is a non-increasing sequence of 
integers placed on this squares. Two partitions p and q are comparable and 
p^qii p[j] ^ q[j] for all j. 

A plane partition P on a Ferrers diagram _F is a non-increasing sequence of 
partitions on the F[j]. The partition on F[j] is noted P[j] and called the j-slice 
of P. Two plane partitions P and Q are comparable and P ^ Q if P[j] ^ Q[j] 
for all j. 

The weight of a partition is its sum; the weight of a plane partition is its sum. 

The plane partitions we are dealing with correspond to the tilings of a 
pseudo-hexagon; they are thus limited, in the sense that the corresponding pile 
must be embedded in the pile that defines the pseudo-hexagon. 
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Definition 4.6 A (s, P) limited plane partition is a plane partition A of the 
integer s such that the plane partitions A and P are comparable and A ^ P. 

We will generate all {s,P) limited plane partitions recursively (see Algo- 



rithm 4.11), adding one j-slice at a time. To this end, we need to generate the 



partitions that are less than a given partition. 

Definition 4.7 A (s, F) limited partition is a partition a of the integer s such 
that a ^ p where p is the partition associated with the Ferrers diagram F. 

Figure ^ (b) shows of a graphical representation of a limited partition. 

The following recursive algorithm uniquely generates all (s, F) limited par- 
tition by adding one by one the parts of the final partitions. The average step 
is thus to add a row with k squares to a yet unfinished Ferrers diagram; the 
value of this part must at least be permitted by the "geometry" of F and be 
less than the last added part. The value must also allow the updated s to fit in 
the remaining "space" ; this is the meaning of r^+i below. 

Algorithm 4.8 

• Input: Integers s,pi, ... ,p„ such that s ^ pi + • ■ • + p„. 

• Output: List C of all (s, (pi, . . . ,Pri)) limited partitions. 

• Compute loop(s,0, (pi, . . . ,Pn)) where 
\oop{t, (ai, . . . ,ag), {pg+i 

— if t = 0, £ ^ (ai, . . 




— else if n = q + 1, C 

— else compute 

loop{t - fc, (ai, ... , a,, fc), (min(pg+2, k),... , min(p„. A;))) 
for k = min(pg_|_i, i, Uq) downto r^+i 
where r^+i = min{l < j < min(pg+i, t, Og) \ 

j + mm{pg+2,j) H \- min(|5„, j) ^ t}. 

Proof of the algorithm A more basic way to write the algorithm would be to 
add any number of squares (fc = l..min(pq+i, t, Uq)) at each step: the min simply 
ensures that the resulting collection of squares is indeed a Ferrers diagram (the 
new part must be less than its counterpart in the Hmiting partition (pg+i), less 
than the number of remaining squares (t) and less than the last added part (a^)) 
and this procedure obviously generates all the desired partitions. 

But it is also wasteful, in that many such combinations would not satisfy t ^ 
Pq+i + ■ ■ ■ + Pn , so we restrict the range of k until its minimal value guarantees, 
by construction, that the following steps will result in a partition limited by p. 
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Since the values we eliminate would not ultimately yield such a partition, we 
still generate all the desired partitions. 

Finally, our algorithm executes a depth-first search on a dynamically gener- 
ated tree T of height n and whose p*'^ level features all the possible combinations 
for the p first rows (the other rows being empty), so that no partition is obtained 
twice. □ 

Complexity Since we perform a depth-first search on T, the execution space of 
the algorithm is bounded by the height of T (n, the number of parts in the limit- 
ing partition) times the number of leaves, which is the number of (s, (pi, . . . ,Pn)) 
limited partitions. We believe no closed formula is known for the latter, but gen- 
erating series techniques (see for instance [ FS96 | and |FS|) yield lemma 4.10| : 



Definition 4.9 The multidegree of a monomial ui*^ •••Mn'" is the sequence 
(ii, . . . , in)- Its cumulated multidegree is the sequence (ii + • ■ ■ + «„, ^2 + ■ ■ • + 



Lemma 4.10 The number of (s,p = (pi,... ,p„)) limited partitions is the 
number of monomials whose cumulated multidegree is less than the conjugate 
partition of p in 

n 

Proof (Sketch) The exponent of Uk in a monomial ui'^ • • • Up'" in P(ui, . . . , Un) 
counts the number of rows of length k in the corresponding partition (see Fig- 
ure |l^). This monomial encodes a partition limited hy p if and only if it satisfies 
ir + ir+i + • • • + *n ^ m(r) -|- m(r -I- 1) -|- • • • -I- m(n) for I ^ r ^ n where m{k) 
is the number of occurences of k in the conjugate partition of p. □ 

Figure 13: Generating series and Ferrers diagrams 



Finally, the execution time of the algorithm is proportionnal to its execution 
space, the multiplying factor being the time needed to compute rg^i at each 
step, which is at most n x pi (exhaustive search). □ 

Note that since truncated generating series can be computed by standard 



programs such as Maple, lemma 4.10 provides another (sub-optimal, but handy) 
way to generate limited partitions. 

We now undertake to generate all the tilings of a pseudo-hexagon, or rather 
all plane partitions limited by a given plane partition P — (Pi, ... , Pn) where 
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the P,'s are the j-shces of P. It is enough to generate all (s,P) Hmited plane 
partitions and to let s range from to the weight u'(P) of P. 

Algorithm 4.11 

• Input: An integer s and a plane partition P = (Pi, . . . ,Pn) such that 
s ^ w{P). 

• Output: The list C of all (s, P) limited plane partitions. 

• Compute loop(s, 0, (Pi, . . . , P„)) where 

l00p(t,(Al,... ,Ag),{Pg+l,... ,P„)) = 

- if t = 0, £ ^ (Ai, ... ,Ag, 0^_^ ); 

(n—q) times 

— else if n = q + 1, 

* generate all (t, min(Ag, Pg+i)) limited partitions with Algorithm 



4.8 



* for each such partition p, C ^ (yli, ... , Aq,p); 
else for all partition k such that 

'w{k) ^ t 
k ^ min(A^,P,+i) 

k + u'(min(fc, Pq+i)) + • • • + ii;(min(fc, P„)) ^ t 

compute 

loop{t-w{k),{Ai,... ,74g,A:), (min(fc,Pg+2), . . . , min(fc, P„))); 
Proof and co mpl exity of the algorithm The analysis is completely analogous 



to Algorithm 4.8's. The execution space is at most n times the number of (s, P) 
limited plane partitions (for which we believe no closed formula is known, but 
as above a characterization with generating series techniques can be obtained) 
and the execution time is at most the execution space times n times the number 
of (u,Pi) limited partitions for u — 0..'w{Pi) (exhaustive search). □ 



Note that Algorithm 4.S is a special case of Algorithm 4.11 , in the same 
manner that a partition is a special case of plane partition; we presented both 
for the sake of clarity. 

5 Domains, fracture lines and seeds 

Up to now, we have considered polygons, which are hole- free (simply connected) 
domains on which one knows how to define a consistent height function, an 
essential ingredient of Section 0. The subsequent results in this paper hold 
provided the results in Section pfdo, but this does not necessarily mean one has 
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(a) A path delimiting a 
domain with a hole 



(b) This domain can be tiled 
even though it has a hole 



Figure 14: Domains with holes 



to deal only with simply connected d omains (see Sectio n 5.1). It would be quite 
exciting to examine how works like Fou97| or | FouOl| extend to the results of 



this paper. 

For the rest of this paper, we consider domains for which the results in 
Section ||hold. To reflect this fact, we will not use the word "polygon". 

A domain T) is defined by its contour path, which is simply any finite-length 
closed path V in the triangular grid. Fracture lines (see Section |5.2D will allow 
us to break the domain into several parts (called fertile zones) whose tilings can 
be generated independently. The lattice of the tilings of V can then be obtained 
as the product of the lattices of the tilings of the fertile zones. Finally, seeds 
(see Section |5.5D will allow us to decompose fertile zones into pseudo-hexagons, 
for which much is already known. 

5.1 Domains with holes 

A domain defined by a finite-length closed path V in the triangular grid need not 
be hole-free: it suffices that V cross itself in a non-trivial way (see Figure |lj (a) 
for an example). Such a domain may still be tileable (see Figure ^ (b)). The 
definitions and results in this paper apply to such domains provided the results 
of Section ^ hold. 

5.2 Fracture hnes, fertile zones and the fracture algorithm 

When considering a general (tileable) domain I?, it may happen that a given 
tile is "fixed", in the sense that it cannot be changed by a flip, so we'd like 



to remove such a tile and concentrate on "interesting" zones (see Section 5.3). 
The tool for this is the fracture line, made of vertices whose height does not 
vary between the minimal and maximal tilings. It allows us to "remove" fixed 
lozenges but there's more: it also allows us to define disjoint sub-domains which 
can be tiled independently (see Section [S^ ), so that the lattice of the tilings 
of T) can be obtained by computing a product of smaller, simpler lattices (see 



Section 5.4). This section attempts to clarify all of the above, step by step, until 
a complete algorithm can be written, proved and analyzed. 

Fracture lines have been investigated in the case of dominoes, notably in 



[Fou96|. In all this section, V is assumed to be a tileable domain such that the 
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results of Section || hold. 

Definition 5.1 Let T be any tiling of V and let h be the associated height 
function. A path P in 2? is T -valid if it is a single vertex of V oi \i\h{v') — h{v) \ = 
1 for u, v' two consecutive vertices of P {i.e. it follows the edges of lozenges 
in T). 

A path P in I? is a T- valid cycle if it is a closed, repetition-free path in T) 
and a T- valid path. 

A subset F of I? is a sub-domain of T) if there exists a tiling T oi T> and a 
T-valid cycle C such that C is the contour path of F. 



Definition 5.2 Let V be any tileable domain and v one of its vertices. We 
denote by Ah(v) the difference between the height of v in the maximal tiling of 
V and its height in the minimal tiling. If Ah{v) = then v is said to be solid. 



Definition 5.3 (Fracture line) A fracture line in 2? is a T-valid cycle in V 
(for some tiling T of V) that is composed only of solid points. Every domain I? 
has at least one such fracture line, namely its contour path, which we will call 
the trivial fracture line of V. 



We n ow undertake to prove Theorem |5.15| . We first prove (see Proposi 



tion 5^) that a solid point in the interior of V cannot be isolated: it belongs at 



least to a solid lozenge. 

Lemma 5.4 Let F be a sub-domain of V, delimited by a valid cycle in a tiling 
T of v. If a triangle t belongs to F, then so does the lozenge £ to which t 
belongs in T. 

Proof Let a, b, c be the vertices of t, and let d be the fourth vertex of £. 
There are exactly two couples made of vertices in {a, 6, c} which have a height 
difference equal to 1 in T. Let a and b be the two vertices which have a height 
difference of 2 (in absolute value), so that [a;b] is not a valid path in T. Since 
F is delimited by a T-valid cycle, [a; b] cannot be part of the contour path of 
F. Therefore d belongs to F, whence the triangle abd belongs to F, and since 
i is the union of the triangles abc and abd, we are done. □ 

Proposition 5.5 Any sub-domain of a tileable domain is itself tileable. 

Proof Let F, a sub-domain of a tileable domain V, be delimited by a valid 
cycle L in a tiling T of 2?. 

If T is a single point, we are done (using no tile). Otherwise, let s be a 
segment of L; it belongs to 2 triangles in the triangular grid, one of which must 



belong to F. By Lemma 5.4, the lozenge £ to which this triangle belongs in T 
belongs to F. If T = £, we are done. 

Otherwise, let Ti be the subset of T) obtained by removing £ from F. It 
is a priori composed of several disjoint subsets of 2?, but each of these subsets 
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is delimited by a T-valid cycle so that the situation is equivalent to a single 
domain Fi, which is then a sub-domain of V. 

Let Fp be the domain obtained by recursively removing a lozenge until noth- 
ing can be done. If Fp contained a triangle it would also contain a lozenge by 



Lemma 5.4 and so it would not be a fixed point. Therefore it contains no trian- 
gle, which means that it is a single point and is a collection of lozenges and 
therefore tileable. □ 

Corollary 5.6 A fracture line of V delimits a tileable sub-domain of T). 
Proof Indeed, a fracture line is a valid cycle so it delimits a sub-domain. □ 

Lemma 5.7 Let si, . . . , sg be six vertices of V that define a hexagon of side 1. 
Let c be the center of this hexagon. If Ah{sk) 7^ for k = 1..6, then Ah{c) 7^ 0. 

In other words, if c is a vertex in the interior of V that cannot be flipped, 
then at least one of its neighbours cannot be flipped. 

Proof Let ft-min and ft-max denote the minimal and maximal height functions. 
Sinc e c and Sk belong to the same triangle, /imin(c) ^ /iinin(sfc) -I- 2 (see Sec- 
L). Moreover /imax(sfc) ^ ^min(sfe) + 3 since Ah{sk) 7^ 0, from which we 



tion 



3.1 



derive /imin(c) ^ /imax(sfe) — 1. In the maximal tiling, no vertex can be a local 



minimum (see Proposition |3.7| and Definition B.8), so ft.max(c) cannot be less 



than /imin(sfc) — 1 for all k — 1..6, whence /imax(c) 7^ /imin(c). □ 

Proposition 5.8 Let V be any tileable domain, not limited to a single point, 
and let w be a solid vertex of V. Then v belongs to a fracture line that encloses 
at least a lozenge. 

Proof The result is trivial if v is on the boundary path of V. Wc now assume 



that we are not in this case. By Lemma 5.7 we know that there exists a solid 
vertex v' at distance 1 (in the triangular grid) of v. Let T be any tiling of 7) 
and let h be the associated height function. If \h{v) — h(v')\ is equal to 2 then 
the lozenge whose middle segment is [v; v'] is defined by two fracture points and 
we are done. 

If \h(v) — h{v')\ = 1 then {v,v') is a valid path in T. Let now si, . . . ,55 
and v' be the (distinct) vertices of the hexagon around v. We prove ab absurdo 
that at least one of the Sk is solid: assume that none of them is solid, so that 

hminisk) ^ /imax(Sfe) for fc = 1..5. 

• Let us assume that h{v) — h{v') — 1. We also have h{v) < /imin('Sfe) + 2 ^ 
(/iinax(s*:)^3)-|-2 ^ /imax (**: ) ~ 1 SO that h{v) would bc less than the height 
of any of its immediate neighbours, i.e. it would be a local minimum, and 
thus Ah{v) could not be 0. 

• Let us now assume h{v) — h{v') + 1. We also have h{v) ^ /imax(sfe) — 2 ^ 
{hmin{sk) + 3) — 2 = ft.min(s*:) + 1 SO h{v) would be morc than the height 
of any of its immediate neighbours and could be down-flipped. 
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Figure 15: The union of two pre-fertile zones need not be pre- fertile 



Either case brings a contradiction so at least one of the Sk is solid; let us 
note it v". If \h{v) — h{v")\ = 2 then, as above, we are done. Otherwise, we have 
shown that v has two immediate neighbours which are solid and linked to v by 
an arc in T. The same result applies to v' and v" so that by induction, since 
the number of vertices in T) is finite, one obtains a fracture line of V. Since the 
domain delimited by this line contains at least a triangle, it contains at least a 



lozenge by Lemma 5.4. □ 



We have now completed our preliminary study of fracture lines. Before 



proving Theorem 5.15, we need to know a bit more about the sub-domains 
containing no fracture lines. There are two cases, according to whether or not 
all the vertices of the sub-domain are on its contour path. We start with the 
second case. 

Definition 5.9 A sub-domain V of D is pre-fertile if: 

• it is connected and tileable; 

• at least one of its vertices is not on its contour path P; 

• Ah{v) 7^ for any such vertex; 

• two vertices of P are at distance 1 in the triangular grid if and only if they 
are neighbours in P. 



A sub-domain is pre-fertile if all its inner points can be flipped at least once. 
The union of two pre-fertile zones need not be a pre-fertile zone. Consider for 
instance Figure |l^ where Vi is the left-hand side hexagon and I?2 is the right- 
hand side one, so that A belongs to the contour path of both. One sees that 
A cannot be involved in any flip, whence Ah{A) — 0; and yet it is not on the 
contour path of Pi U I?2- For this reason we give the following definition (both 
Vi and I?2 are examples of such zones) : 

Definition 5.10 (Fertile zone) Let I?i and 2?2 be two pre-fertile sub-domains 
of T>. Their union is fertile if Vi U I?2 is itself pre-fertile. A sub-domain of D is 
fertile if it is pre-fertile and a maximal element for fertile union. 



Proposition 5.11 Let F be a pre-fertile sub-domain of T). It is fertile if and 
only if its contour path is a fracture line of V. 
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Proof We first prove that the contour path of a fertile sub-domain _F is a 
fracture line in V. Let a be any vertex of the contour path of F. If a belongs to 



the contour path of P, it is solid by Lemma 2.7. Otherwise, a is in the interior 
of P; let Ha be the hexagon (in the triangular grid) of side 1 whose center is a. 
Since F is fertile it is a maximal element for fertile union so that F U Ha is not 
pre-fertile, which means that at least one of its inner vertices is solid. Since 
none of F's interior is, a must be this vertex. In other words, all the vertices on 
the contour path of F are solid and we are done. 

Conversely, we prove that a pre-fertile sub-domain whose contour path is a 
fracture line in V is fertile. Let F' be any other pre-fertile sub-domain of V 
such that F n F' — 9 and U F' is connected. Let a be a vertex that belongs 
to both F and F' (it is on the boundary path of both). If a is not an inner 
point of U F' for all such a, then this union is not pre-fertile (see the last part 



of Definition 5^). Otherwise, a is in the interior of U F', and it is also solid 
since it belongs to the boundary path of F, whence F U F' is not pre-fertile, 
from which we deduce that F is indeed a maximal element for fertile union and 
therefore a fertile sub-domain. □ 

Proposition 5.12 Let F be a sub-domain of D, delimited by a fracture line L, 
such that at least one of its vertices is not on its contour path. Then F is fertile 
if and only if it contains no non-trivial fracture line of V other than L. 

Proof If F is fertile then Ah{v) ^ for any of its inner vertices, so that it 
can contain no fracture line other than L. Conversely, recall that a fracture line 
can be made of only one point: therefore if F contains no fracture line other 
than L th en A h{v) for any of its inner vertices. Moreover F is tileable (see 
Corollary |5.6| ), connected, and two vertices of L that are not neighbours in L 
cannot be at distance 1 in the triangular grid (this would define a non-trivial 
fracture line) so it is pre-fertile, and since its contour path is a fracture line of 2? 



it is a fertile zone by Proposition 5.11, □ 



We have completed our study of fertile zones; we now turn to sub-domains 
whose vertices are all on the contour path. 

Lemma 5.13 Let L be a fracture line of T) enclosing a domain F, not limited 
to a single point, such that all the vertices of F belong to L. Then either F is 
a lozenge or it contains a fracture line which defines a lozenge. 

Proof Since F is not limited to a single point, it contains at least a triangle. 
This triangle belongs to a lozenge a; in a tiling T of V. The fracture line L is a 



T- valid path, therefore x is embedded in F by lemma 5.4. By hypothesis all its 
vertices (vi, V2, V3 and 114) are on a fracture line of V, so that Ah{vk) = for 
fc = 1..4. Thus the closed path i;i-ti2-W3-W4-'yi is a fracture line in V and we are 
done. □ 

We need one last definition, which corresponds to the sub-domains that 
are obtained by recursively looking for fracture lines, and therefore used by 
Algorithm 5.16 (see also Figures |l| (b) and0 (b)): 
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Definition 5.14 A fracture zone in a tileable domain D is a sub-domain of "D 
not reduced to a single point, delimited by a fracture line and containing no 
other fracture line. 



Theorem 5.15 (Fracture theorem) Let D be any tileable domain such that 
the results of Section ^ hold. A fracture zone in V is either a fertile zone or a 
lozenge. 

Proof Let _F be a fracture zone of V. If all its vertices belong to its contour 
path, then it is a lozenge by Lemma 5.13| . Otherwise Ah{v) ^ for any inner 
vertex v of F by Pr oposition 5.8 since _F is a fracture zone, so that F is fertile 
by Proposition |5.12| . □ 

We now have the tools to extract the fertile zones from a domain: 
Algorithm 5.16 (Fracture algorithm) 

• Input: A tileable domain T). 

• Output: The fertile zones of V. 



• Step 1: Use Thurston's Algorithm 2.£ to build the maximal and minimal 
tilings of v. Store the solid points of I? in a list S. 

• Step 2: Use S and the minimal tiling of V to build all the fracture lines 
except the trivial one; store them in a list L. 

• Step 3: Remove from L all the fracture lines consisting of exactly 4 
vertices. 

• Step 4: Return L. 



Proof of the algorithm Everything has already been proved in Theorem 5.15 



Complexity Steps 1 and 2, which control the execution time of the algorithm, 
are both linear in the number of vertices of V. 



5.3 An example of the use of the Fracture Theorem 

We now present an example. Let us consider the finite- length closed path of 
Figure [l^ (a), which delimits a tileable domain (see Figure (b) for its minimal 
and maximal tilings). Using fracture lines, one can disconnect the central part, 
whose tilings give rise to a trivial lattice (see Figure |l^ (a)). Thus only the 
fertile zones remain (see Figure 0(b)). 
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(a) A generic (b) The minimal and maximal 

domain tilings of the domain 

Figure 16: A domain and its extreme tilings 





(a) The maximal triv- (b) The fertile zones in the minimal 

ial fracture zone and maximal tilings 

Figure 17: The fracture zones 



5.4 An application of the Fracture Theorem 

Since fertile zones are delimited by fracture lines, flips done inside one of them 
cannot have any effect on the other fertile zones: therefore one can study in- 
dependently the tilings of these zones. Since the tilings of V will be obtained 
by "gluing" tilings from these disjoint zones, their lattice can be obtained by 
product. Consider Figure in which the domain contains two fertile zones. 
We will see later (see Section |^) an algorithm to generate the lattice of the 
tilings of such a zone; on this example each fertile zone is a pseudo-hexagon. 



so its tilings can be found using Algorithm 4.11 . We build the product lattice 
of these lattices by connecting every element in the left-hand side lattice with 
every element in the right-hand side lattice, which finally gives the lattice of the 
tilings of the initial domain. 

5.5 Seeds and their fiUings 



We have seen in Section 5.2 that fertile zones play a major role in the generation 



of all the tilings of a generic domain V. We have also seen in Section 5.4 that 
we can restrict our study to these zones and later compute a product lattice. 
We now show that fertile zones can be decomposed into collections of pseudo- 
hexagons (see Figure |l^ (a) for a start). The correspondence between tilings 
and compact piles of cubes allows us to use geometric terms; in other words, an 
up-flip can be viewed as adding a cube. 

Definition 5.17 (Seed) A seed is the minimal tiling of a hexagon of side 1. 
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The minimal 
tiling of a do- 
main 



The lattices of 
the tilings of the 
fertile zones 



How to compute 
a product lattice 




The complete lattice 
of the tilings of the 
domain 



Figure 18: The main steps in computing the lattice 



Note that a pseudo-hexagon (see Definition |4.3| ) contains exactly one seed. 
We now investigate the immediate properties of seeds. 



Lemma 5.18 The union of two pseudo- hexagons sharing the same seed is again 
a pseudo-hexagon. 



Proof Using Definition 4.3, it is enough to show that the union of two compact 
piles of cubes, aligned on the same set of axes, is itself a compact pile of cubes. 
But this immediately follows from Definition 4.1. □ 





(a) A range and the maximal range (b) Covering a seed with a cube po- 
of a seed tentially creates three new seeds 



Figure 19: Seeds and ranges 



Definition 5.19 The range r{s) of a seed s in a tiling T of I? is the union of 
the pseudo-hexagons that contain it. The maximum range R{s) of s is the union 
of its ranges when T ranges among all the tilings of V. 

The set of pseudo-hexagons surrounding a seed s is not empty since it 
contains at least s itself. Thus the range of s is well-defined according to 



Lemma 5.18 
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Definition 5.20 A tiling of r(s) (resp. the maximal, minimal tiling of R{s)) 
will be called a filling of s (resp. the maximal filling, minimal filling of s). 
We denote by Max(s) and Min(s) the maximal and minimal fillings of s; by 
extension, Max(I?) and Min(I?) are the maximal and minimal tilings of V. 

The intuitive idea behind the term "filling" is that one gradually fills r(s) 
with (compact piles of) cubes. Note that a filling of a seed is a compact pile 
of cubes and that the maximal filling of s is obtained by doing all the possible 
up-flips in r{s). 

Adding a cube to cover a seed s (in other words, flipping s) generally creates 
three new seeds si, S2 and S3 as shown in Figure ^ (b). We thus see that 
there is a natural partial order on the ranges of seeds, defined by inclusion. In 
our example, r(s,) C r(s) for i G {1,2,3}. We will be mostly interested in the 
maximal elements for this partial order and now proceed to define this more 
precisely. 

Definition 5.21 A seed s is a child of a seed s' if R{s) is embedded in R{s'). 
A seed is a proper seed if it is the child of only itself. 

Proposition 5.22 Every cube in a tiling of V belongs to a filling of a proper 
seed. 

Proof Every cube is the maximal tiling of a hexagon of side 1 and it is added 
when the minimal tiling of this hexagon is flipped. But this minimal tiling is a 
seed, whose range is included in the range of a proper seed. Therefore the cube 
belongs to a filling of this proper seed. □ 

Proposition 5.23 The maximal fillings of two distinct proper seeds are dis- 
joint. 

Proof Let C be a cube belonging to the maximal filling of two proper seeds 
si and S2- Since a filling of a seed is a compact pile of cubes, there is an 
uninterrupted line of cubes connecting C to the base planes defined by si and 
S2- Since these planes vary two by two by a translation, the two basis are in 
fact the same, so that si = S2- D 

Corollary 5.24 The maximal fillings of the proper seeds in a domain 2? form 
a partition of the subset of Z"^ defined by the maximal tiling of the fertile zones 
oiV. 

We know that the set of the tilings of T) is connected by flips; therefore 
maximally filling all the seeds that appear in Min(I?) potentially creates new 
flippable zones: that is, new seeds and new proper seeds (see Figure ^). To 
clarify the situation, we give the following definition: 

Definition 5.25 A proper seed is of order if it appears in the minimal tiling 
of T>. It is of order n + 1 if it is a proper seed in the tiling of V obtained when 
all the proper seeds of order fee {1, • . . , are maximally filled with cubes. 
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Two seeds of One seed of Two lozenges around The neighbourhood 
order order 1 a seed of order 1 of a seed of order 1 

Figure 20: Seeds 



5.6 Seeds of order 1 

Let us examine the basic properties of a seed s of order 1. Remember that this 
is the minimum tihng of a hexagon and that it does not appear in the minimal 
tiling of v. Thus at most two of the three lozenges that make up s can come 
from Min(I?). 

Let us first assume that exactly two of the three lozenges in s come from 
Min(I?), and exactly one from a proper seed t of order 0. Then s appears while 
filling the range of t, and it is a flippable zone: therefore it is not a proper seed 
and should be flipped while filling t; it is not a seed of order 1. 

So at least two of the lozenges in s come from fillings of proper seeds of order 
0. For reasons of symmetry one can assume that they are as in Figure Since 
they come from fillings they belong to cubes, so that the immediate neighbour- 
hood of s looks like the right-most picture in Figure |2^. 

5.7 Seeds of greater order 

For a given domain V the order of any seed is bounded since Min(I?) and 
Max(X') vary by a given number of cubes and each proper seed contributes at 
least one cube. Conversely, if n is any integer, there exist domains in which one 
can find proper seeds of order n. To illustrate this, consider Figure |2l|. If one 
maximally fills all the seeds of order (b), a sub-domain closely resembling the 
first can be outlined (c); indeed, only the size differs. 

Thus one can build a family (I'm) of domains in which covering all the seeds 
of order of exactly yields I?„ and therefore all the seeds which are of 

order k in are of order fc + 1 in I?ri+i- Using this procedure, one can generate 
seeds of any order. 

6 C-minimal tilings and intervals in jC{V) 

A close inspection of the lattice C{T>) of the tilings of T> will reveal that it 
contains interesting intervals, which we will use to generate all the tilings of a 
domain by performing the required fiips in an orderly fashion. 
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(a) Seeds of or- (b) Seeds of or- (c) Seeds of order 1 in I?„ 

der in I?„ der 1 in I?„ are of order in 

Figure 21: A family of domains 






(a) A lozenge in Max(r') and the cor- (b) Min((7,)(T) with #{C,} = 2 
responding C-minimal tiling 



Figure 22: Tilings defined by the presence of a lozenge 



6.1 C-minimal tilings of a pseudo-hexagon 

Let V he a pseudo-hexagon and let C be a cube in a filling of the proper seed s 



of I?; as we will see in Section 3.S, marking such a cube can be very useful. 
What can one say about the fillings of s in which C appears? Since fillings are 
compact piles of cubes, each such filling must also contain all the cubes in the 
parallelepiped defined by the diagonal going from s to C (see Figure ^ (a)). 

Definition 6.1 Let P be a pseudo-hexagon and C a cube in a filling of its 
proper seed s. The C-minimal tiling of P (which we denote by Minc(-P)) is 
the tiling associated with the filling of s obtained by starting from Min(T) and 
adding all the cubes in the parallelepiped defined by the diagonal (s, C). 

This definition is consistent. First, a C-minimal tiling is a tiling since the 
parallelepiped {V{C) ) is a compact pile of cubes. Second, if a filling of s contains 
7'(C) then it contains C. Finally, if a filling of s contains C then it must contain 
7'(C) since a filling is a compact pile. 

Definition 6.2 Let P be a pseudo-hexagon and (Ci)i^i^„ a collection of cubes 
in a filling of the proper seed s of P. Let Tk be the Cfc-minimal tiling of P for 
k € {1, . . . ,n}. The {Ci)-minimal tiling of T (which we denote by Min((7.)(P)) 
is Sup(Ti, ... ,r„). 

An example of such a tiling is given in Figure |2^ (b) . 

All of the above can be reformulated in a dual way. Instead of adding cubes 
from the minimal tiling, remove cubes from the maximal one: this can be viewed 
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Figure 23: The successive V'' for k — 0..d{'D) 



©3 



as adding "anti-cubes^\ the rule being that cubes and anti-cubes annihilate each 
other. We could thus define a C-maximal tiling as a tiling of I? which does not 
contain the cube C . The notion can be extended to families, thus defining 
[Ci]- maximal tilings. 

6.2 The fundamental intervals of CiV) 

Let now T) be any domain (such that the results of Section ^ hold) and let T{T>) 
be the set of its tilings and C{T>) the associated lattice. We have seen that 
Max(I?), seen as a pile of cubes, can be partitionned according to the maximal 
fillings of its proper seeds. 

Definition 6.3 We denote by 2?", n ^ 1, the tiling of P obtained by maximally 
filling the ranges of all the proper seeds of order fc g {1, . . . , n}. By convention, 
T)'^ = Min(2?). The degree d{T>) of T) is the minimum of the integers n such 
that = Max(I?) . A cube C is of order k if it belongs to the filling of a seed 
of order k. 

Note that Max(P'=) Min(X''=+i) for ^ fc < d{V). See Figure H for an 
example. 

Notation 6.4 The set of all the tihngs t of V for which =<; t =^ P'^+i is 
an interval in the lattice of the tilings of V; we denote this interval by Tn{T>), 
^ n < d{T>)). By extension, T(X>) denotes the set of all the tilings of T). 

Definition 6.5 The 7^(1?), ^ n < d{'D), are the fundamental intervals of the 
lattice C{p) of the tilings of T). 

Proposition 6.6 (7o(I?),... , 7^(u)_i(I?)) is a maximal chain of intervals in 
C{V). 

Proof It follows from the definitions that = Min(X'), X''^^^) = Max(X') 
and Maxp'^) = Min[V^+^) for fc < d(V). □ 

A graphical representation of this chain is given in Figure |^. 
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-T{V) 



Ti{V) 



%{V) 



Min(I') 



Figure 24: The fundamental intervals of T(P) 





(a) The contour path of a C- 
minimal tiling is a range of 
a proper seed 





(b) Building the smallest of the 
tilings containing C 



Figure 25: C-minimal tilings in T) 



6.3 C-minimal tilings of a fertile zone 

Let I? be a domain such that the results of Section || hold. Using fracture lines 
(see Section 5^), we can suppose that I? is a fertile zone. Let C be a cube 
that appears in a filling of a proper seed s of order n of T). It is clear from the 
definitions that C G I?". C induces a C-minimal tiling of i?(s), let us denote it 
by P{C), whose contour path is a pseudo- hexagon (whose proper seed is s) and 
therefore a range of s (see Figure p5| (a)). 

Definition 6.7 The C-minimal tiling of V, let us denote it by Minc(r'), is the 
infimum of all the tilings of V that contain C. If (Ci)i^i^„ is a family of cubes, 
the (Ci)-minimal tihng of V is Min(Ci)(2?) = Sup(Minci (V), ... , Minc„ {V)). 

We now sketch a recursive construction of Mine (I?) (see Figure ^ (b)). 
7'(C) marks cubes that belong to and therefore to the maximal fillings 

of seeds ti of order k ^ n — 1. These cubes give rise to (Cj)-minimal tilings in 
the ti, which in turn mark cubes in I?""^, and so on. Since we always choose 
minimal tilings in the ranges of the proper seeds, the tiling is less than any tiling 
containing C, and it contains C, therefore it is Minc>(I?). 



32 



7 An algorithm to generate the lattice of the 
tilings of a general domain 

In Section |4.2| we have proposed an algorithm to generate efficiently the tilings 
of a pseudo-hexagon. We now proceed to the general case, where no hypothesis 
is made on the shape of the domain V, except that the results of Section ^ hold. 

In order to make use of already found tilings and reduce the computation 
time, we need the following important fact: if (C;) is a collection a cubes, all 
of them of order n, the interval of between Min((7.)(I?) and T>" plus the 

cubes Ci is isomorphic to the interval between Min((7.)(D) minus the cubes Ci 
and V". This is obvious because cubes of order k < n can be added to 
Min((7.)(I?) whether the cubes of order n are present or not. 

Formally stated, there is an isomorphism between the intervals [Min((;7.) {V) ; 
Sup(Min(c.)(2?),2?")] and [Inf(Min(c,) (2?), I?") ; 2?"] oi C{V). 



Algorithm 7.1 

• Input: A finite-length closed path V, enclosing a domain V, in the trian- 
gular grid. 

• Output: The lattice C{T>) of the tilings of T>. 



• Step 1: Use Algorithm 5.16 to break the domain into mutually indepen- 



dant fracture zones (see Definition 5.14). The lattice of the tilings of a 



single lozenge is trivial. For each of the fertile zones, follow steps 2 to 4. 

Step 2 {V is now a fertile zone): Recursively build the tilings V'', ^ 
k s; diV). Set the fist L-r, of the tilings of V to 0. 



• Step 3: For fc = to dcg(X') - 1 do: 

to 



— Step 3.1: For each proper seed of order fc, use Algorithm 4.11 
generate all the fillings of the seed. 

— Step 3.2: Compute the cartesian product of these fillings. 

— Step 3.3: Each element of this product is a family (Ci) of cubes; for 
each of these families, compute Min(c.) (2?). 

— Step 3.4: Removing the cubes Ci from Minj^p.-) {V) yields a collection 
of cubes of order at most fc — 1, which we denote by Mm'^^\^{V), or 
the empty set if fc = 0. Look up in Lx> the (already found) tilings 
between Mm^~^^{'D) and I?*^; for each of these tilings t, add t + (Ci) 
to Lx). 



• Step 4: For each tiling T in L-p, compute the height function and find the 
local minima: they correspond to cubes that can be added. Connect T to 
each of the tilings obtained from it by adding exactly one of these cubes. 
This generates the lattice of the tilings of the fertile zone D. 
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step 5: Compute the product of the lattices found for each fracture zone. 



An practical implementation of this algorithm should update the lattice of 
the fertile zone at each newly found tiling. 

Proof of the algorithm/. Lot T be any tiling of a fertile zone. It can be associated 
with a unique pile of cubes. Let Af be the maximum of the orders of these cubes. 
The cubes of order M belong to the fillings of the seeds of order M , and all of 
them have been taken care of by construction. The cubes of order at most M —1 
in T form a pile of cubes which has been generated in step 3.4. We conclude 
that T has been encountered at least once by the algorithm. Moreover, all the 
tilings produced by the algorithm are distinct. Therefore we have generated 
exactly once each tiling of each fertile zone of V, and hence exactly once each 
tiling of P. □ 

Space complexity: Let 1^(1?) | denote the number of tilings of the domain D (we 
believe no closed formula is known). Each tiling of 2? is generated only once; 
the number of links in £(2?) starting from a particular vertex is at most the 
width of the lattice of which we believe nothing is known except that 

it is trivially bounded by |T(2?)|; and there is a small overhead to account for 
the V'^. The execution space of the algorithm is thus 0(|T('D)p). □ 

Time complexity: The execution time of the algorithm is controlled by steps 3 
and 4, but since the average (or worst case) number of seeds and size of their 
maximal ranges is (totally unknown and) highly dependent on the shape of the 
domain, no non-trivial bound can be given for the time being. A rough analysis 
can be conducted as follows. In step 3.3, one needs to compute a minimal 
tiling for a familiy of cubes, of which there is at most the number of cubes in 
the maximal tiling of V, which in turn is less than the number T(T>) of tilings 
of V (since each cube is bijectively related to the associated minimal tiling). 
The minimal tiling for each cube can be computed in a number of steps that is 
again at most the number of cubes in the maximal tiling of V. Looking up the 
already found tilings in step 3.4 can be done in at most |T(I?)| operations. Since 
steps 3.3 and 3.4 must be conducted for each tiling, the overall cost of step 3 
is 0{\T{T>)\^). In step 4, computing the height function requires I'D] operations 
for each tiling and connecting a tiling to its fathers in the lattice requires at 
most T(V) operations, so that the global cost of step 4 is 0(|T(I?)p x |X>|). 

Thus an upper bound for the time complexity is 0(|T(D)|^), but the execu- 
tion time is probably far less in practice. □ 



8 Conclusion and perspectives 

We have shown that the isomorphism between Conway and Lagarias' lozenge 
group and Z'^ is quite fruitful; it justifies the intuitive geometric interpretation 
of lozenge tilings and allows us notably to define the proper seeds and the 
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important intermediate tilings T)^ ^ of which Thurston's minimal and maximal 
tilings of V are particular cases. 

Also, we have made much use of the lattice structure; this seems to indicate 
that lattice theory is a tool well adapted to the study of tilings. 

We believe that the techniques developped in this paper may be extended 
to study some related questions: 

• Enumerating the tilings of D, without generating them, is a difficult ques- 



tion; Algorithm 7.1 can be used to exhibit an exact (but not very useful) 
formula, which is turn could provide good bounds if one were able to count 
the number of plane partitions that are both greater and smaller than two 
given plane partitions. Cruder bounds can be obtained using the lattice 
structure. 

In the case of tilings with dominoes there is a natural definition of flips, 
and therefore of seeds. The range of a seed s in a tiling T can be defined 
as follows: mark all the seeds in T except s as forbidden (they should not 
be flipped) and perform all the possible up- flips. The range of s is the 
collection of squares around the vertices that have undergone a flip in the 
process. It seems therefore likely that our algorithms could be adapted 
straightforwardly to the case of dominoes. The geometric interpretation 
might be preserved using Levitov tiles (see | RYOOf ) . 



Our definition of (the maximal tiling of) a pseudo-hexagon is a natural 
generalization of a Ferrers diagram since it is merely a geometrical rep- 
resentation of a plane partition. Let us now attribute an integer to each 
cube in such a way that sequences are non- increasing along each axis: this 
is a generalization of a Ferrers diagram in dimension 4 (this is called a 



solid partition in [Macl6|). Thus one can easily define a generalization 



of pseudo-hexagons in dimension p ^ 2 and may therefore be able to 
study tilings in this dimension. In particular, the recursive algorithms of 



Section i.1 should be readily upgradable. 

If more counting results were known, one could use Algorithms 7.1 and 4.11 
to generate tilings uniformly at random. 

The C-minimal tilings seem to correspond to the sup-irreducible elements 
of the lattice of the tilings. It would therefore be interesting to study 
the link (with Birkhoff's theorem) between these tilings and the order 
associated with the lattice. 
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